Galilean Equivalence for Galactic Dark Matter 
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Satellite galaxies are tidally disrupted as they orbit the Milky Way. If dark matter (DM) experi- 
ences a stronger self-attraction than baryons, stars will preferentially gain rather than lose energy 
during tidal disruption leading to an enhancement in the trailing compared to the leading tidal 
stream. The Sgr dwarf galaxy is seen to have roughly equal streams, challenging models in which 
DM and baryons accelerate differently by more than 10%. Future observations and a better under- 
standing of DM distribution should allow detection of equivalence violation at the percent level. 
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Galileo showed at the Leaning Tower that objects of 
different masses and materials fall the same way in a grav- 
itational field. This equivalence principle, later a corner- 
stone of Einstein's general relativity, has been tested re- 
peatedly through a variety of experiments since Galileo's 
time and since the time of Einstein (see, e.g., Ref. 0). 
These tests have confirmed that the materials found on 
Earth and in the Solar System all satisfy the equivalence 
principle to a remarkable degree. But what about the 
dark matter that fills galactic halos and dominates the 
mass density of the Universe? Is the dark matter in the 
Milky Way's halo accelerated the same as baryons in a 
gravitational field? 

The simplest and most favored candidates for dark 
matter (DM), like weakly interacting massive particles 
[S El and axions H, EI , do satisfy the equivalence 
principle (EP). However, there are a number of reasons to 
test this assumption. First of all, we still have no empiri- 
cal evidence for the existence of WIMPs or axions. Some 
have argued that a stronger self-gravity for dark matter is 
required to clear dwarf galaxies from voids in the galaxy 
distribution 0, 0, ^j. Moreover, the recent discovery 
that the cosmological expansion is accelerating 0, 
suggests that there may be more to gravity than general 
relativity — in particular, the quintessence field may me- 
diate an additional long-range self-interaction between 
dark-matter particles (l3l,ll4j . There is thus considerable 
motivation to scrutinize our cherished notions about the 
equivalence principle. 

Violations of the equivalence principle in the dark sec- 
tor may be modeled phenomenologically by attributing 
to dark-matter particles ip a "fifth force" [lM Hij , 



9 



Anr 



(1) 



Here g is a dimensionless coupling constant, and is the 
mass of the scalar particle <fi mediating the interaction. 
On scales r <C m^ 1 , the potential of Eq. leads to an 
inverse-square-law force between DM particles of mass 
m^p with a strength suppressed by a factor /? 2 compared 
to gravity, where f3 = gmpi/y/^nm^. 



Several cosmological consequences of such a DM force 
have already been explored. First of all, to clear dwarf 
galaxies from voids, values /3 > 1 and mT 1 > 1 Mpc are 
required |l7j . An attractive force for r -C m^ 1 would 
enhance structure formation on these scales leading to a 
corresponding increase in the density-perturbation power 
spectrum |l8j . an effect, though, that can be mimicked by 
a blue tilt in the power spectrum. An EP-violating cou- 
pling between DM and quintessence could also induce a 
scale-independent bias between baryons and DM,_though 
this effect is model dependent ^4|. Refs. fl5t fl8| had 
noted that a DM force would strip a baryonic core from 
its dark halo, and applied this to typical galaxies in the 
Coma cluster to set a limit < 2.2. Clusters might also 
test an attractive DM force, as baryons would be pref- 
erentially lost compared to the more tightly bound DM 
during the mergers leading to their formation [19(. This 
test is complicated, however, by gas physics which is ex- 
pected to reduce the cluster baryon-to-DM mass ratio 
below the cosmological value, even in the absence of a 
DM force. 

In this Letter, we consider the effects of a DM force 
on galactic scales. We propose here that tidal streams 
produced by the disruption of a DM-dominated satellite 
galaxy orbiting in the halo of a much larger host galaxy 
provide a powerful probe of an EP-violating DM force. 
The reasoning follows by comparing the satellite's orbital 
energy E OI b, the energy £" t id imparted during tidal dis- 
ruption, and the self-binding energy Ebm of the satellite 
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Here the host galaxy has a potential <&host(-R) and mass 
Mr within the satellite's orbit of radius R, and the satel- 
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lite has a mass m sat and radius r sa t which fill its tidal 
radius ryd- When the satellite is much less massive than 
the host galaxy, m sat /M R Cl, a distinct hierarchy, 
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exists in these three energy scales, implying that the dis- 
rupted stars and satellite will trace similar orbits in the 
host galaxy's potential regardless of the details of tidal 
disruption or the satellite's internal structure. The dis- 
rupted stars will act like purely baryonic test particles, 
while the satellite itself behaves largely like a DM test 
particle, if it is DM dominated. 

Fortunately, the Sagittarius (Sgr) dwarf galaxy, the 
Milky Way's closest satellite at a Galactocentric dis- 
tance of only 16 kpc, is nearly ideal for our purposes. 
The Sgr dwarf has extended leading and trailing tidal 
streams observed by the Two-Micron All Sky Survey 
[2MASS) [HI and the Sloan Digital Sky Survey (SDSS) 
|22j. Using a sample of over 1,000 M-giant stars with a 
known color-magnitude relation, the 2MASS collabora- 
tion have measured not just surface brightnesses along 
the streams, but distances and spectroscopic velocities 
as well [ijjj. Comparing these observations to simula- 
tions has led to estimates of the mass of the Sgr dwarf of 
M Sgr = (2-5) x 10 8 M Q , mass-to-light ratio M Sgr /i Sgr = 
14 — 36 Mq/Lq, and Sgr orbit with pericenter 10-19 kpc, 
apocenter 56-59 kpc, and period 0.85-0.87 Gyr [24,]. The 
large mass-to-light ratio suggests that the Sgr dwarf is 
indeed DM dominated and therefore a suitable place to 
search for DM forces. 

To study more carefully the effects of EP violation on 
tidal disruption, we performed our own simulations of the 
tidal disruption of a satellite with a mass (5 x 10 8 M Q ), 
mass-to-light ratio (40M Q / Lq), and orbit (pericenter 14 
kpc, apocenter 59 kpc) similar to that of the Sgr dwarf. 
We could not compare our simulations directly with those 
of Ref. ji^l, as we performed TV-body simulations of a 
Navarro-Frenk- White (NFW) profile for our Milky Way 
halos, and they used a static logarithmic potential. An 
active halo allows for dynamical friction over the course of 
the simulation and possible backreaction on the halo due 
to the DM force. While we did not attempt to reproduce 
the detailed features of the Sgr tidal streams, our simula- 
tions are sufficient to demonstrate that even a small DM 
force could have significant observational consequences. 
The initial conditions for our simulations were produced 
using GALACTICS 25], which makes use of phase-space 
distribution functions (DFs) that are analytic in the or- 
bital energy and angular momentum. By Jeans' theo- 
rem, these DFs are equilibrium solutions to the collision- 
less Boltzmann equations [i^ , and they can be combined 
to produce realistic and stable models of the composite 
Milky Way bulge-disk- halo system pffi ] . We used the two 
Milky Way models of Ref. |25j that best fit observational 
constraints, including the Galactic rotation curve and lo- 
cal velocity ellipsoid. The simulations were evolved using 
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FIG. 1: Simulations of the tidal disruption of a satellite galaxy 
in the presence of a dark-matter force. The charge-to-mass 
ratio j3 increases from 0.0 in increments of 0.1 going counter- 
clockwise from the bottom left. The Galactic disk is in black. 
Sgr stars are shown in red (dark grey) while the Sgr dark mat- 
ter is blue (light grey). The tidal streams are projected onto 
the orbital plane. Orbits are counterclockwise; the upper left 
figure shows that for f3 = 0.3 (a dark-matter force 9% the 
strength of gravity) stars are almost absent from the leading 
stream (at 12 o'clock with respect to the Galactic center). X's 
denote the location of the bound Sgr core. 



a modified version of the TV-body code GADGET-2 [2jJ. 
A more detailed description of our simulations are pro- 
vided in Ref. [H. 

Four simulations of tidal disruption are depicted in 
Fig. H with DM forces given by Eq. 0) with different 
values of the charge-to-mass ratio (3. The scalar field is 
assumed massless (m^ = 0), so the DM force is a true 
inverse square law. The ratio (3 increases from 0.0 at 
bottom left to 0.3 at top left as one proceeds counter- 
clockwise. The simulations begin with the satellite at 
apocenter 59 kpc from the Galactic center and last for 
2.4 Gyr (almost three full orbits). The tangential veloc- 
ities are adjusted so that all orbits are projected to have 
a pericenter of 14 kpc. The orbits are counterclockwise 
in the x-z plane, so that the edge of the leading stream 
appears at 12 o'clock with respect to the Galactic center 
in Fig.^ while the edge of the trailing stream is at about 
10 o'clock. The Sgr dwarf is modelled with a truncated 
NFW profile for both stars and DM, in keeping with the 
simulations of Ref. [3] , where it was concluded that ob- 
servations could not yet determine distinct profiles for 
the two components. Thus, the stars shown in red (dark 
grey) in the bottom left panel are simply a downsampling 
of the DM distribution illustrated in blue (light grey). 

As the DM force increases in strength, the leading 
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FIG. 2: Surface density of stars as a function of angular dis- 
tance 9 along the tidal stream. The satellite core is located 
at 0°, while the trailing and leading streams are at positive 
and negative respectively. As the tidal streams wrap around 
the Galaxy more than once, 9 extends beyond ±180°. The 
four curves correspond to the four panels of Fig. with 
black (solid), red (long-dashed), green (short-dashed), and 
blue (dotted) curves belonging to the f3 = 0.0 through 0.3 
simulations. 



stream is systematically depleted of stars, while the trail- 
ing stream is correspondingly enhanced. By the time (3 
reaches 0.3 in the top left panel, the leading stream is 
virtually devoid of stars. The primary reason for this 
effect is that in the presence of an attractive DM force, 
the center of mass of the satellite's stars is displaced out- 
wards with respect to that of its DM. The bound stars 
lie at the bottom of the satellite's gravitational potential 
well and are therefore forced to orbit the Galactic center 
at the same speed as the DM. However, they do not have 
the attractive pull of the DM force from the Milky Way's 
halo to supplement gravity in providing the required cen- 
tripetal force. The stars are therefore displaced outwards 
so that the inward gravitational pull of the satellite's DM 
can provide this additional centripetal force. From this 
outer position, stars are more likely to be tidally dis- 
rupted from the far side of the satellite than the side 
closest to the Galactic center. Stars disrupted from the 
far side gain energy and are boosted onto higher orbits 
in the Milky Way's potential well where their angular 
velocity is slower than that of the satellite. They there- 
fore trail behind the satellite and develop into a trailing 
tidal stream. A repulsive DM force will induce an op- 
posite effect, displacing the stars towards the Galactic 
center and preferentially creating a leading rather than a 
trailing tidal stream. 



This asymmetry in the leading compared to the trail- 
ing tidal streams is a distictive signature of a DM force 
that can be observed in the stellar densities measured 
along the stream. The normalized stellar densities of 
the four simulations presented in Fig. [J are shown in the 
four curves of Fig. [3 Orbits in our composite Milky Way 
model do not close, and the four peaks in the stellar den- 
sity distribution correspond to the four apocenter pas- 
sages appearing at 2, 7, 10, and 12 o'clock in Fig. ^ Or- 
bital velocities are minimized at apocenter so stars tend 
to accumulate there. The ratio of the number of stars 
near the apocenters furthest along the leading and trail- 
ing tidal streams thus provides a convenient measure of 
the asymmetry between the streams. The ratio of the 
number of stars in the leading segment stretching from 
—300° to —200° as compared to the trailing segment from 
350° to 450° drops from 0.66 in the absence of a DM force 
down to 0.44, 0.091, and 0.0042 for = 0.1, 0.2, and 0.3 
as indicated by the solid black curve in Fig.01 The SDSS 
has observed hundreds of stars per square degree along 
the Sgr tidal stream [22| ■ As the Sgr dwarf is observed to 
have an extensive leading stream, we conclude that a DM 
force as weak as 9% the strength of gravity is likely to be 
observationally unacceptable, in which case the proposal 
EH E3 that EP-violating dark matter clears dwarf 
galaxies from voids would be ruled out. 

Our simulations suggest that current observations can 
already place impressive constraints on a DM- force, but 
several concerns remain before we can confront our simu- 
lations with data. If the stream wraps around the Galaxy 
more than once, we must be able to distinguish true lead- 
ing stars from trailing stars that have almost been lapped 
by the satellite. We have been able to do this surpris- 
ingly well in simulations using only the radial velocities, 
distances, and positions along the stream. As 2MASS 
has collected this data, identifying leading and trailing 
stars should already be feasible and can ceratainly be 
accomplished by a future high-precision astrometry ex- 
periment like the Space Interferometry Mission (SIM) or 
Gaia. More troublesome is whether some other change 
in our Sgr or Milky Way models could produce the same 
asymmetric tidal tails that we are claiming as a signa- 
ture of a DM force. Future investigation of this concern 
is certainly needed, and we have made a first attempt at 
this in Ref . |28| , the results of which are summarized in 
Fig. |3 This signature is seen to be robust to changes in 
the Milky Way model and the mass, orbit, and phase- 
space distribution of the Sgr dwarf. In the absence of a 
DM force, the leading-to-trailing ratio exceeds 0.5 for all 
our models, while for > 0.2 the ratio is always below 
0.2. The detailed morphology of the stream also allows us 
to anticipate when the ratio should be high or low, mak- 
ing our test more sensitive than a single number. Our 
simulations suggest that it may be possible to detect a 
DM force a few percent the strength of gravity. We may 
not be able to drop DM off the Leaning Tower of Pisa, 
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FIG. 3: The ratio of leading to trailing stars as a function 
of charge-to-mass ratio j3 for different models of the host- 
satellite system. The black (solid) curve is our default best- 
fit model. The magenta (long-short dashed) curve doubles 
the initial mass of the satellite. The two red (long-dashed) 
curves have rotating satellites: the top is prograde and the 
bottom is retrograde. The blue (dotted) curves have satellites 
with different orbits: the top curve has a more circular orbit, 
while the bottom curve has a planar orbit rather than the 
polar orbit of Sgr. The cyan (dot-short dashed) curve uses 
a Milky Way model with lighter halo and heavier disk. The 
green (short-dashed) curve has a satellite where 25% of the 
most bound particles represent stars. The yellow (dot-long 
dashed) curves have satellites with lower M/L ratios (higher 
stellar mass fractions): the top curve has M/L = 4.5 while 
the bottom curve has M/L = 10. 

but the Sgr tidal streams may be the next best thing. 
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